DOC  FILE. COPY, 


NRL  Report  (237 


iD  ' 

CJ> 

<£> 


Equations  of  Powered  Rocket  Ascent 
and  Orbit  Trajectory 


M.  H.  Reilly 

Communication  Sciences  Division 


1 furs^ 

a l,«  my  t-3  a y 


May  7,  1979 


D D C 

icrsitr-m 

JUN  4 1319 

A 


NAVAL  RESEARCH  LABORATORY 
WnUagim,  D.C 

Approved  for  public  release;  distribution  unlimited. 


an  fun '•••' 


■*?»  i.n-jufWEUWf 


^CLAtyriCATiQN  or  this  raok  figw  pw« 

1 REPORT  DOCUMENTATION  PAGE 


„ II.  OOVT  ACCESSION  NO. I 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
UujiCNT7!  CATALOQ  NuSilN 


s.  "title  iw  jStttn) 


y „ EQUATIONS  OE,  POWERED  ROCKET  ASCENT  AtfD 
-1]  "ORBIT  TRAJECTORY  4 / 


I rvFS  UF  MFO««  * FCNIOO  COVCNCO 

/ Final  repSkt.dui  one  phase  of  a 
^continuing  <4ltL  problem 
a.  performing  one.  report  nmmskm 


ir.  AuTftO"t«; 


(/*J>  M:H.|R^UIy 


r fcnfoamino  onoanhation  name  and  aod.tcss 

Naval  Research  Laboratory 
Washington,  D.C.  20375 

■■  ■ " «— « ■ ' ■«  ' 1 1 < w ■ 

it.  controlling  opficc  name  ano  adorcsi  / 

Communication  Sciences  Division  ( / / 

Telecommunications  Environmental  Effects  Branch 


I.  CONTNACT  on  ONANT  NUMlCNrij 


PNOOAAM  element,  fnojcct,  task 
Aar  a a won*  unit  auMacas 


7^K03-7T 


MONlToAlNO  AOEnCV  NAME  a AODftESSfN  dlllerent  leetyjCimttrelllrtg  Olllco)  I*.  SCCuaiTV  CLASS,  (o I Hilt  teporl) 

f ) •)  Y 7/V^..^  UNCLASSIFIED 

\ ' y/  ' ' ‘ , If*.  OCC^  ASSiriCATICN/oOWNOaADINO 


If.  DlSTaiaUTlON  STATEMENT  (el  idle  Report) 


Approved  for  public  release;  distrihutipn  unlimited. 


[ IT.  OlSTNieuTiON  STATCMCNT  (el  the  ed tired  metered  In  Bleed  19.  tl  dHlptent  Item  Report) 


If.  SlinFLEMENt  ANY  NOTES 


[ 19.  KEY  IHOROI  fCwil/nu*  on  r»Fnr«i  if  n«e«««*.Y  md  Identity  by  Woe*  nubiarl 


Rocket  trajectory 
Dynamic* 

Earth’s  curvature 
Rotation  effects 


Orbit  injection 

Satellite  orbits  and  trajectory 


10.  ABSTRACT  fCcnffmi*  on  r«v«ro  »ld*  it  nm c**»*fy  And  Identify  by  b/oefc  numborj 

^Equations  of  rocket  and  satellite  dynamics  are  obtai-  ed  at  various  levels  of  approximation 
and  make  possible  analyses  of  powered  rocket  ascent,  orb  t injection,  and  subsequent  orbital 
motion.  Rocket  and  satellite  trajectory  are  specified  by  these  equations,  which  include  the  effects 
of  earth's  rotation  and  curvature.  Various  analysis  techniques  are  given,  so  that  the  investigator 
can  choose  the  simplest  one  suited  to  his  needs.  In  the  process,  material  not  readily  available  is 
presented,  and  some  problems  of  ambiguity  and  inaccuracy  in  standard  references  arc  resolved,  iv 


00  1 JAn*7J  1473  coition  or  I NOV  at  IS  OMOLCTC 


S/N  0 10 }-0 14- 6601 


tfCUNITV  CLASSIFICATION  O f THIS  FAOC  (Mim  Dote  Inteted) 


1 gsi  9SO 


lr  1,L 


mm**'*  B;uf;r.--»:i;<.-rT-  -»y  1 


CONTENTS 


v 


r ' 7*  vV.  ".*  tcnZ _ 

» — i 1 .-  bi  ’ i.t.V  CotlRS 


Avail  and/or 

Mat 

k 

special 

1.  INTRODUCTION  1 

2.  POWERED  ROCKET  ASCENT  NEGLECTING  EARTH’S 

ROTATION  EFFECTS  2 

‘The  Equations  of  Motion 2 

Integration  of  the  Equations  of  Motion  6 

3.  INCLUSION  OF  EARTH’S  ROTATION , 7 

4.  ORBIT  TRAJECTORY  FROM  INJECTION  CONDITIONS 

AT  ROCKET  BURNOUT 16 

5.  VEHICLE  FLIGHT  SEEN  BY  INERTIAL  AND 

EARTH  OBSERVERS 20 

6.  DISCUSSION 21 

7.  ACKNOWLEDGMENTS  22 


8.  REFERENCES  22 

APPENDIX  A - Lagrangian  Method,  Including  the  Presence  of 

Nonconservative  Forces  23 


APPENDIX  B - Derivation  of  the  Lagrangian  and  Generalized 

Forces  for  the  Rocket  Problem 25 

APPENDIX  C — Velocity  Increment  Estimates  for  Various  Effects 

in  First-Stage  Rocket  Motion  27 

APPENDIX  D - Determination  of  Launch  Profile  Parameters  from 

Desired  Orbit  Injection  Conditions  in  a I lat  Earth  Approximation 28 

APPENDIX  E - Thrust  Increase  and  Aerodynamic  Effects  30 


EQUATIONS  OF  POWERED  ROCKET  ASCENT 
AND  ORBIT  TRAJECTORY 


1.  INTRODUCTION 

In  the  analysis  of  systems  in  space,  one  often  encounters  the  problem  of  specification  and 
design  of  a launch  vehicle  and  its  trajectory,  so  that  certain  performance  requirements  may  be 
met.  Also  of  interest  are  the  flight  path  of  the  space  vehicle  in  orbit,  particularly  with  respect 
to  interactive  ground  systems,  and  the  relation  of  the  orbit  to  the  launch  profile.  As  will  be 
explained  below,  however,  it  is  usually  difficult  to  extract  the  requisite  information  from  vari- 
ous standard  references.  This  report  should  help  to  remedy  this  situation.  One  of  the  primary 
goals  of  this  report  is  to  present  a logical  exposition  of  approximations  and  sufficient  informa- 
tion and  guidance  to  facilitate  an  investigator's  choice  of  the  simplest  analysis  technique  suited 
to  his  needs.  He  may  thus  possibly  be  able  to  avoid  an  unnecessarily  complex,  time-consuming 
analysis  technique,  such  as  the  "full”  numerical  solution  approach  to  the  problem. 

In  many  standard  reference  sources  on  powered  rocket  ascent  from  the  earth,  equations 
are  introduced  in  an  ad  hoc  and  incomplete  manner.  There  are  also  inaccuracies  and  incon- 
sistencies to  contend  with.  For  example.  Ball  and  Osborne’s  Eqs.  (1-20)  and  (1-22)  (1]  leave 
out  certain  kinematic  terms  which  are  related  to  a radial  gravitational  field  and  could  become 
important  for  high  rocket  velocities.  Ehricke  (2]  correctly  introduces  the  effect  of  these  terms  in 
his  Eq.  (5-27)  in  an  ad  hoc  fashion,  but  then  he  seemingly  incorrectly  introduces  their  effect 
later  in  his  Eqs.  (6-39)  and  (6-40),  also  in  an  ad  hoc  fashion.  Ruppe  [3]  gives  a basically 
correct  treatment  of  these  terms  in  his  Eqs.  (3.1M3.4),  but  from  the  outset  his  treatment 
ignores  the  orientation  or  angle  of  attack  of  the  vehicle  through  his  use  of  a gravity  tilt  condi- 
tion. Furthermore,  Ball  and  Osborne  [1]  give  a sketchy  treatment  of  earth's  rotation  effects,  in 
which  an  initial  rocket  velocity  imparted  by  earth’s  rotation  is  included  in  a flat  earth  approxi- 
mation. Ruppe  [3]  does  a more  complete  job  on  this,  but  his  results  are  also  restricted  to  a flat 
earth.  In  most  instances  a powered  rocket  ascent  will  not  cover  a large  enough  ground  range  to 
necessitate  taking  curvature  of  the  earth  into  account,  but  in  some  cases  it  will.  To  cover  these 
cases,  where  ground  range  might  be  in  excess  of  500  n.mi.,  a correct  treatment  of  earth’s  cur- 
vature effects  is  necessary,  but  it  is  absent  from  the  references  that  have  been  discussed  [1-3]. 

This  report  is  intended  is  to  remedy  the  above  deficiencies  and  present  in  one  place,  in  a 
consistent  notation  set,  and  in  reasonably  coherent  fashion  the  basic  equations  of  powered  rock- 
et ascent  and  orbit  trajectory  needed  for  analysis  techniques  at  various  levels  of  approximation. 
To  start  with,  in  the  first  part  of  Section  2,  we  obtain  the  rocket  equations  of  motion  in  an  iner- 
tial frame.  As  in  nearly  all  elementary  treatments,  the  complications  from  the  fact  that  the  at- 
mosphere is  not  stationary  in  an  inertial  reference  frame  are  initially  ignored.  In  the  iBtter  part 
of  this  section  these  equations  of  motion  are  integrated  to  give  velocity  and  coordinate  expres- 
sions as  a function  of  time  in  an  inertial  frame.  In  Section  3,  we  consider  the  complications  of 
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earth's  rotation,  including  the  effect  of  a corotsting  atmosphere.  A "full*  solution  of  the  rocket 
trajectory  would  thus  involve  a complicated  three-dimensional  numerical  solution  of  the  equa- 
tions of  motion.  In  Section  3 it  is  shown  that  the  great  complications  which  arise  from  this 
course  of  action  can  be  avoided  by  solving  simplified  equations  of  motion  in  the  earth-fixed, 
rotating  reference  frame  for  the  period  of  powered  rocket  ascent  when  aerodynamic  effects  are 
important;  the  data  are  then  transformed  back  to  the  inertial  frame,  which  amounts  to  addition 
of  a rotation  velocity  vector;  and  finally  the  inertial  velocities  computed  as  in  Section  2 for  por- 
tions of  the  rocket  flight  when  aerodynamic  effects  are  not  important  are  added.  The  results  of 
Section  3 demonstrate  that  approximate  velocity  and  coordinate  expressions  for  rocket  flight  can 
be  obtained  by  a rather  simple  extension  of  results  from  Section  2.  The  expressions  in  Section 
3 will  be  found  to  include,  to  a good  approximation,  not  only  earth  and  atmosphere  rotation 
effects,  but  also  earth’s  curvature  effects  for  the  totJ  powered  trajectory.  In  Section  4 the 
nature  of  the  elliptical  satellite  orbit  and  its  point  of  entry,  as  determined  from  the  injection 
conditions  which  exist  at  cutoff  of  the  rocket  engines,  are  considered.  In  Section  5 the  pro- 
cedures are  described  for  determining  the  vehicle  flight  coordinates  versus  time  in  both 
powered  and  orbital  phases,  and  as  seen  by  both  inertial  and  earth  observers.  In  Section  6 there 
is  a discussion  of  the  equations  obtained  and  their  underlying  validity. 

2.  POWERED  ROCKET  ASCENT  NEGLECTING  EARTH'S  ROTATION  EFFECTS 

The  Equations  of  Motion 

The  coordinate  system  used  to  describe  rocket  motion  is  shown  in  Fig.  1 . The  forces  on 
the  rocket  are  the  thrust  force  of  the  engine,  gravity,  and  aerodynamic  forces.  As  in  most  ele- 
mentary treatments  11-3],  we  assume  initially  that  the  atmosphere  is  stationary,  so  that  aero- 
dynamic forces  on  the  rocket  arise  purely  because  of  the  rocket  motion.  Actually,  the  atmo- 
sphere corotates  with  the  earth  to  a first  approximation,  but  for  targe  launch  vehicles  this  turns 
out  to  have  consequences  of  only  secondary  importance  for  rocket  motion  (3).  We  do,  how- 
ever, include  this  as  one  of  the  rotating  earth  effects  in  Section  3.  With  the  assumptiom  of  a 
stationary  atmosphere  and  an  axially  symmetric  rocket  in  our  inertial  frame,  the  simplest  case 
of  rocket  motion  will  be  planar  and  conveniently  describable  in  circular  coordinates  (see.  Fij. 
1),  because  of  the  radial  nature  of  the  gravitational  field  and  the  spherical  earth.  The  effects  of 
earth's  rotation  will  be  disregarded  for  now,  but  will  be  included  in  Section  3.  Shown  in  Fig.  1 
are  the  radius  vector  r from  the  center  of  the  earth  to  the  satellite,  its  associated  altitude  y,  the 
angular  displacement  from  launch  4,  which  gives  the  ground  range  x when  the  radius  of  the 
earth  R is  factored  in,  the  flight  path  angle  of  the  center  of  mass  of  the  rocket  ^ with  respect  to 
the  local  horizontal  (or  "heading"),  and  the  angle  of  attack  a of  the  rocket  axis  with  respect  to 
the  flight  direction. 

For  a - 0 in  Fig.  1 the  only  aerodynamic  force  on  the  symmetric  rocket  will  be  a drag 
force  opposite  to  the  flight  direction.  If  one  defines  drag  D and  lift  L aerodynamic  forces  for 
a * 0 as  being  antiparallel  and  perpendicular  to  the  flight  direction,  it  is  found  for  small  a that 
L is  linearly  proportional  to  a and  that  the  correction  to  D is  approximately  quadratic  in  a [2]. 
Unlike  Ehricke  (2)  and  Ruppe  {3],  Ball  and  Osborne  (ll  unconventionally  refer  drag  and  lift 
direction  to  the  rocket  axis  of  symmetry  (or  "roll*  axis).  In  Fig.  2 wc  indicate  a force  diagram 
for  the  rocket  and  include  a brief  description  of  the  symbols.  We  use  the  standard  assumptions 
found  in  many  references  11-3J: 
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Fig.  I — Coordinate  syslem  for  two-dimensional  rocket  motion 


• - CENTER  OF  GRAVITY  l 

O - CENTER  OF  PRESSURE  W 

w - FLIGHT  PATH  HEAOING  f, 

« - ANGLE  OF  ATTACK 

T - THRUST  FORCE  /j 

/»  - ANGLE  OF  THRUST 

0 - DRAG 


- LIFT 

- WEIGHT 

- VECTOR  FROM  CENTER  OF  GRAVITY 
TO  CENTER  OF  COMBUSTION 

- VECTOR  FROM  CENTER  OF  GRAVITY 
TO  CENTER  OF  PRESSURE 


Fig.  2 — Force  diagram  for  the  rocket  with  symbol  legend 


• the  aerodynamic  forces  act  through  the  center  of  pressure, 

I • the  force  of  gravity  acts  through  the  center  of  gravity,  and 

i 

{ • the  thrust  force  is  applied  through  the  "center  of  combustion". 

| 

• One  may  utilize  the  Lagrangian  method  in  this  problem,  suitably  generalized  to  include 

f the  presence  of  nonconservative  forces  (4).  A brief  review  of  this  approach  is  included  in 

Appendix  A.  In  Appendix  B this  method  is  applied  to  the  present  problem  in  finding  the 
f Lagrangian  L and  generalized  nonconservative  forces. 
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U) 

(3) 

(3) 

(4) 

In  effect.  Bell  end  Otborne  111  leave  out  the  term  rd1  in  Eq.  (2)  end  e term  r d in  Eq. 
O)  in  their  Eqs.  (1*20)  end  (1-22)  (Note  that  r,  rk  correspond  to  their  y.  x).  Equation*  (2) 
end  (3)  ere  the  equation*  of  motion  for  the  ret’  at  and  angular  coordinates  of  the  center  of 
mess,  and  Eq.  (4)  applies  to  the  coordinate  a for  the  internal  motion  (cf.  Fig.  1).  The  center  of 
mess  Eqs.,  (2)  end  (3),  in  which  we  are  particularly  interested,  can  be  manipulated  as  follows. 

We  have  for  the  velocity  of  the  center  of  mass: 

rd  - v cos  d . r - v a in  d.  (5) 


Lagrange's  equations  are  found  as 

± U JIL  , 

where  g , - r,  d.  end  »-•  in  turn.  Their  evaluation  yields 

r - rd*  d K/i3  - {L/M)  cos  d - {DIM)  sin  d d {TIM)  sin  (d^ot  ft). 
rd  d 2rd  - - {L/M)  sin  d - {DIM)  cos  d d ( TIM ) coe  (d  d o d ft),  and 
/a  • IjU  cos  a d D sin  al  - l\  T sin  ft . 


so  that 

v1  - P d r1*1. 


Differentiating  this  relationship  and  using  equations  *.2),  (3),  and  (S),  we 

vv  - rr  4-  r*dd  d rdlr 

- r (rd1  - K/r1  d (L/W)  cos  d - (D/,V)  sin  d d {TIM)  sin 
d rd  (-  2 rd  - (L/M)  sin  d - iD/M)  cos  d d (T/M)  cos  (d  d 

- v sin  dl-  K/r1  d (i/V)  cos  d - (0/M)  sin  d d (7*/M)  sir 
d v cos  d l~  iL/M)  sin  d ~ {D/M)  cos  d d {T/M)  cos  (d 


obtain 

(d  d a d 8)  | 
a d 0)1  d rd*r 
(d  d o + /J)) 
dad  fl)}. 


Combining  terms,  we  have 

v - - (K/r>)  sin  d - (D/M)  d {T/M)  cos  (a  d 3).  (6) 


This  can  also  be  written  as 

v - {T/M)  cos  ft  cos  a ~ {D/M)  - {K/r1)  sin  d - {T/M)  sin  ft  sin  a,  (6') 


and  this  is  the  form  of  Ehricke's  Eq.  (5-23)  12).  To  obtain  the  equation  for  d.  wo  differentiate 
Eq.  (5): 


Hence, 


v cos  d - v (sin  d)d  • rd  d rd- 


- ~ sin  d - •jj  + ^ cos  (a  d ft) 


v (sin  d)d 


cos  d - rd 
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4-  2rd  4-  ■—  sin  * ♦ ■—  cos  * - jj  cos  (4  + o 4-  R). 

Combining  term*  and  using  Eq.  (S),  we  have 

vd  « - - ~|  cos  d 4.  sin  (o  4-  01  ♦ H) 

This  pm  be  written  at 

vd  • ~ coe  0 tin  a 4-  + “ tin  0 coe  a - — ~-|  cos  d,  (7’) 

which  corresponds  to  the  form  of  Ehricke's  Eq.  (5-27)  (2).  The  derivation  of  Eqa.  (6)  and  (7) 
it  very  similar  to  Ruppe’s  procedure  for  obtaining  his  Eqs.  (5.3)  and  (3.4)  (3). 

The  flrst  level  of  simplification  which  is  analyzed  elsewhere  (2.3).  results  from  the  fact 
that  for  the  large  launch  vehicles  of  interest  here.  0 « a.  so  that  the  equations  of  motion 
simplify  to 

v — t/C/r1)  sin  d ~ (D/M)  4-  ( T/M ) cos  a and  (8) 

vp  =»  - [(A/r1)  - (vVr)]  cos  d 4-  (L/M)  4-  ( T/M)  sin  o.  (9) 


Another  great  simplification  results  when  one  finds  (2,3]  that  the  increase  of  thrust  at  higher 
altitudes  is  approximately  cancelled  by  the  drag  effect  in  the  calculation  of  cutoff  conditions. 
Small  a is  normally  a requirement  when  aerodynamic  forces  ore  a factor  (e.g.,  for  first-stage 
motion  or  altitudes  less  than  60  km).  Under  these  conditions  the  lift  can  be  written  as  in  Ref. 
2: 

UQS  pvV2  =*  (dQ/d«)  a S pvV2,  (10) 


where  S is  a reference  cross-sectional  area  of  the  vehicle,  p is  the  atmospheric  density,  and 
BCL/da  is  a constant  determined  from  wind-tunnel  testing.  Hence,  at  this  level  of 
simplification, 

v =«  - (Af/r1*  sin  4 + (T/M)  cos  « (11) 


and 


COS  d + 


r 

M 


sin  a, 


(12) 


one  can  do  relatively  simple  calculations.  Small  a ascent  is  desirable  from  the  standpoint  of 
launch  efficiency,  i.e.,  minimum  fuel  expenditure  for  accomplishing  the  mission  (2,3],  Thrust 
is  considered  to  be  a constant  in  Eqs.  (11)  and  (12),  with 

r-Aifv,.  (13) 
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Here  v,  is  the  effective  speed  with  which  exhaust  gases  are  ejected  relative  to  the  nozzle  exit 
aperture  (.1),  and  M is  the  rate  at  which  exhaust  mass  is  ejected.  One  simplified,  approximate 
procedure  for  equations  (10  and  021  is  to:  (a)  assume  a programmed  deflection  function  +01 
for  the  various  launch  vehicle  stages,  which  results  in  correct  orbit  entry  conditions  at  burnout, 
(b)  soUe  Eq  0 0 by  integration,  and  (cl  compute  angle  of  attack  from  021.  The  last  step  in 
this  procedure  is  simply  to  check  the  assumption  of  reasonably  small  n used  in  obtaining  Eqs. 
(Ill  and  021  In  this  context  « is  given  as  the  solution  of 


Mv»  + M (K/r}  1 U - vrVA)  cos  » 
T (1  + (dC',/dol  S^/2  r) 


041 


This  equation  turns  out  to  be  identical  to  the  result  derived  by  Ehricke  (cf.  his  Eqs.  (6*401,  (6* 
421,  and  (b-SOll,  exert*  for  the  square* bracketed  term  in  Eq,  (141,  which  Ehricke  refers  to  as  a 
centrifugal  load  factor  (2).  He  introduces  this  factor  in  an  ad  hoc  fashion,  and  he  apparently 
incorrectly  replaces  v in  this  factor  by  v cos  •.  Ehricke  shows  that  for  a large  winged  rocket 
vehicle  the  term  enclosed  in  braces  in  the  denominator  substantially  reduces  angle  of  attack  a 
beneath  its  vacuum  value,  particularly  for  the  first  stage  of  the  launch  vehicle 


Normally,  a is  small  in  Arst-stage  motion  when  aerodynamic  effects  are  important,  but 
can  be  fairly  large  in  higher  stage  motion  when  aerodynamic  forces  are  insignificant,  depending 
on  required  orbit  injection  conditions.  The  presence  of  the  cos  a term  in  Eq.  (Ill  can  thus  be 
important  for  higher  stages. 


Integration  ef  the  Equations  ef  Maths 

We  follow  Ehricke's  notation  |2]  in  integrating  Eq.  (Ill  for  a particular  stage  of  a multi- 
stage rocket  vehicle.  We  introduce  r,  as  the  burn  duration  of  the  st.  je  and  the  normalized 
time  variable  ( for  that  stage: 

i - f /!,.  (151 

Here  f varies  between  0 and  1.  During  this  time,  mass  varies  as 

M/M,  - (Ma  ~ Mt)/M,  - 1 - ~ - 1 - (A  , (161 

"■0 

where 

A 3 MtJM,  • W,/W,  (171 

is  s parameter  which  is  given  as  the  weight  of  propellant  for  that  stage  divided  by  the  initial 
weight  of  the  launch  vehicle  at  that  stage.  If  we  assume  small  angles  of  attack,  such  that  cos 
« = 1,  and  normalize  the  vehicle  speed  to  the  exit  speed  of  the  exhaust  ga»es,  i.e.,  define 

X s v/v,.  (181 

wo  find  er.ily  (cf.  Eq.  <6-24>  in  Ref.  1211  that 

Av/v,  - Ax  u x ~ X,  - -Ml  - <Al  - (A/it.1  (/,!«.  (191 

where  n,  :.s  the  thrust-to-weight  ratio 

n,  - T;M,giR)  - TfW,.  (201 

where  g (rl  - K/r*  is  the  acceleration  of  gravity,  and 
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(/.)<  « f*  singd*  (21  > 

Mere  4 varies  from  its  value  at  the  beginning  of  the  particular  stage,  which  is  the  same  as  its 
cutoff  value  at  the  end  of  the  preceding  stage  (intermediate  coasting  can  be  included  as  a 
separate  stage),  to  its  value  at  the  particular  time  according  to  a well-defined  prescription  for 
d(<)  Ehncke  uses  a deflection  program  such  that  d vanishes  at  ( - 0 and  I and  ♦ - 0 at 
( - I (21  This  is  not  necessarily  optimum.  The  integration  in  Fq.  (21)  is  performed  by  a sim- 
ple numerical  method,  and  so  *({)  in  Eq  (19)  is  determined  at  the  particular  ( mesh  points 
selected  for  the  numerical  integration.  Ruppe  |3)  gives  a procedure  for  obtaining  analytic 
renults.  Then  from  Fig.  1 and  Eq.  (S),  altitude  and  ground  range  are  determined  from 

y ~ V,  - v,r,  x «n 

* - x.  - v,r,  J*  x cos  * - Rfd  - d.)  (22) 


From  Eq.  (14).  the  angle  of  at  taw  is  found  from 

WHIN* 

ana*  


cos  4 +■*-(!-  i\)  Q 

•V  df 


do  IT 


(23) 


The  initial  parameters  , v and  y,  for  the  particular  stage  in  question  are  the  same  as  the 
cutoff  values  for  the  preceding  stage.  For  the  very  first  stage  - >•„  • 0,  and  we  would  use 
X»  • 0 in  our  calculations.  We  should  reiterate,  however,  that  earth's  rotation  effects  have  not 
been  included.  Fortunately,  the  additional  complications  caused  by  their  inclusion  do  not 
present  unsurmountahle  problems,  as  we  will  attempt  to  demonstrate  in  the  next  section 

3.  INCLUSION  OF  EARTH'S  ROTATION 

The  preceding  equations  would  be  valid  in  an  inertiel  frame  of  reference  if  all  effects  of  a 
rotating  earth  could  be  excluded.  One  effect  of  a rotating  earth  is  to  impart  an  initial  velocity  to 
the  rocket  at  launch.  If  this  were  the  only  effect,  this  velocity  could  be  simply  vectorially  added 
to  the  velocities  found  in  the  preceding  section  and  integrated  to  give  a correction  to  altitude 
and  ground  coordinates.  There  is  another  effect,  however;  the  earth's  atmosphere  corotates 
with  the  earth  to  a first  approximation  (disregarding  the  ordinary  winds  experienced  by  an  earth 
observer).  If  the  atmosphere  were  stationary  in  an  inertial  reference  frame,  an  earth  observer 
would  he  subjected  to  constant  wind  speed  of  903  cos  L knots,  where  L is  the  terrestrial  latitude 
of  the  observer.  This  suggests  that  an  exact  solution  of  the  rocket  problem  might  profitably  be 
carried  out  in  « frame  of  reference  which  rotates  with  the  earth,  at  least  while  aerodynamic 
effects  are  important.  This  ts  the  basts  of  an  approximation  suggested  by  Ruppe  (3).  For  the 
sake  of  this  discussion,  we  suppose  that  the  first  stage  of  our  multistage  rocket  vehicle 
corresponds  to  the  period  of  rocket  flight  when  aerodynamic  effects  are  impottant.  Ruppe  sug- 
gests {31  that  the  equations  of  motion  should  be  solved  in  the  rotating  reference  frame  (in 
which  the  atmosphere  is  stationary)  for  the  first  stage,  that  the  results  should  be  transformed 
back  to  the  inertial  frame  of  reference,  and  subsequent  stages  u .oed  in  the  inertial  frame, 
since  aerodynamic  effects  are  unimportant  for  the  higher  stages.  V«»e  transformation  from  the 
rotating  reference  frame  back  to  the  inertial  reference  frame  after  cutoff  of  the  first  stage  gives 
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At,  - At*  | + t.,  (2*) 

where  At,  it  the  rocket  velocity  in  the  inertial  frame  after  cutoff  of  stale  1.  and  Av*,  it  the 
tamo  velocity  determined  in  the  rotating  reference  frame.  From  the  general  theory  of  tuch 
transformations  (tee,  e.|.,  Kef.  S) 

*#*l!i*fi*  «r*|l  ♦ "jjrj  cot  Lt&  a v,£.  (25) 

where  r,  it  the  radius  vector  to  the  satellite,  and  y,  its  altitude  after  cutoff  of  the  Aral  stale. 
The  approximation  indicated  in  Eq.  (25)  is  associated  with  a "flat  earth*  approximation  for  the 
Arst-ttage  motion  in  which  the  curvature  of  the  earth  is  disregarded,  to  that  L,  is  the  latitude  of 
the  launch  point,  and  t is  a unit  vector  in  the  easterly  direction  at  the  time  and  place  of  launch. 
This  approximation,  which  is  discussed  later,  is  very  often  valid  for  Arst-stage  motion,  and  for 
the  modest  altitudes  attained,  one  sees  that  v,  is  approximately  the  velocity  imparted  to  the 
rocket  at  its  launch  point  by  earth's  rotation  in  an  inertial  frame.  The  magnitude  of  mcR  is  t5] 
4o5  m/t  (1 524  ft/s)  “ 903  knots.  The  true  inertial  velocity  during  the  Ath  stage  is  given  by 

t(r,  ♦ r , f ...  + fi4)  - v,  4>  At g,  + Atj  + • • • + (Av4)c  * t,  + v'.  (26) 

Here  Avy0  > D it  the  velocity  increment  of  the  / th  stage  (computed  u in  the  previous  sec- 
tion), and  (A>k)(  it  the  velocity  increment  in  the  *th  stage.  The  entity  ry  is  the  burn  duration 
time  for  the  j th  stage,  and  ( is  the  reduced  time  variable  which  varies  between  0 and  1 (cf.  Eq. 
(15)).  We  shall  And  Atft,  later  and  justify  that  it  and  v'  in  Eq.  (34)  are  approximately  coplanar, 
just  as  in  the  previous  section.  But  for  now,  we  simply  assume  this  fact  and  proceed  lu  carry 
out  the  vector  addition  of  t0  and  v’  in  Eq  (26).  We  do  this  with  the  help  of  Fig.  3.  which 
shows  the  planar  motion  represented  bv  v and  the  vector  t„  drawn  from  the  launch  point  in  an 
easterly  direction.  The  orientation  of  the  plane  of  v’  is  specifted  by  the  azimuth  angle  a.  meas- 
ured from  north  (N)  to  east  (E)  at  the  launch  point.  The  vector  v'  is  speciAed  by  its  magnitude 
v'  and  by  its  direction,  which  can  be  deternpned  from  a,,  the  heading  A',  and  the  ground  range 
angle  A’.  all  of  which  are  indicated  in  Fig.  3.  To  facilitate  the  vector  addition,  we  break  v0  into 
a component  vaJL  perpendicular  to  the  plane  of  v'  and  a component  v,,, , in  this  plane.  Evi- 
dently, 


V#1  - va  cos  a,  and  v*, , - v,  sin  a*.  (27) 

Hence,  t can  be  speciAed  as  the  sum  of  three  perpendicular  component  vectors  (shown  in  Fig. 

4), 


t - vr  + v*  + v#a.  (28) 

where  v,  is  parallel  to  the  radius  vector  from  the  center  of  the  earth  to  the  rocket,  and  v*  is 
perpendicular  to  it,  but  in  the  plane  of  v'.  From  the  geometry  in  Fig.  4 it  follows  that 

V*  - v’  cos  A'  + v„  |COS  A’. 

v,  - v’  sin  A‘  + v„ , sin  A'.  (29) 

and 


V - (v*  + V*  + 

The  true  heading  A with  respect  to  the  local  horizontal  is  given  by 

sin  A " v„/v.  (30) 
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EQUATOR 


Fig  3 — Vector  orientations  for  inclusion  of  earth's  rotation  effects 
in  a Pat  earth  approximation  for  first  stage  (see  test) 
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Fi«.4  — Vector  diagram  for  inclusion  of  earth's  rotttion  in  a 
flat  earth  approximation  for  first  stage 

In  order  to  find  the  azimuth  angle  a for  the  vehicle  flight  path  projection,  we  need  to  make  use 

of  spherical  trigonometry  formulae.  If  we  disregard  the  effect  of  initial  velocity  from  earth’s 

rotation,  we  can  relate  the  entities  in  Fig.  S by  the  use  of  the  laws  of  cosines  and  sines: 

sin  L'  - sin  f*  cos  4‘  + cos  L,  cos  a0  sin  4' 

sin  a'  - sin  a,  cos  L,/ cos  L' 

sin(/'  - 4)  ■ sin  a,  sin  $’/cos  L'.  (31) 

These  account  for  the  variations  of  latitude,  azimuth  angle,  and  longitude,  respectively,  in  the 
plane  of  v’  [31.  Now  from  Fig.  6, 

v*  sin  a'  + vol  cos  a'  v*  tan  a'  + voI 

tan  q ■■  j ■ "j  ■ * — ■ ■ r . (32) 

v h cos  a - voX  sin  a v*  - v0l  tan  a 
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Fig.  6 — Vector  diagram  Tor  azimuth  angle  determination 
taking  into  account  earth’s  rotation 


We  thus  have  completed  the  vector  addition  of  Eq.  (26).  With  the  knowledge  of  v.  and  a, 
we  can  determine  the  actual  position  of  the  space  vehicle  in  an  inertial  coordinate  system.  We 
can  describe  these  position  coordinates  as  the  altitude  of  the  space  vehicle  plus  the  latitude  and 
longitude  of  its  flight  projection  on  the  surface  of  the  earth  in  an  inertial  geocentric-equatorial 
coordinate  system  [S].  (In  this  system  the  latitude-longitude  grid  is  fixed  in  space.)  With 

X s v/v„  (33) 

we  have 

y -y„  - v,r,  * sin  t </f 

1 - Lo  - vi  //  *£ ff cos « w 
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I - 4 - v,/,  /;  ®in  a rff  (34) 

for  the  changes  in  altitude,  latitude,  and  longitude  during  a particular  bum  stage  of  the  rocket 
vehicle.  The  parameters  in  the  above  integrals  are  evaluated  for  values  of  ('  between  0 and  £, 
where  f — 0 is  associated  with  cutoff  of  the  preceding  stage,  and  ( - 1 is  associated  with 
cutoff  of  the  present  stage. 

The  preceding  expressions,  particularly  those  which  invoke  the  approximation  in  Eq.  (25), 
can  be  considered  valid  in  a flat  earth  approximation  for  the  first  stage— the  period  of  rocket 
ascent  when  aerodynamic  effects  are  significant.  This  approximation  for  any  part  of  rocket 
motion  can  be  considered  valid,  for  example,  when  4'  can  be  disregarded,  as  in  Eqs.  (29)  and 
(31).  We  note,  for  example,  that  a ground  range  of  500  n.mi.  corresponds  to  =*  8.3*.  For 
this  value  sin  4'  is  only  about  1/7  the  value  of  cos  4'  - 0.99,  so  that  sin  - 0 and  cos  - 
1 is  an  adequate  approximation,  particularly  for  the  significant  portion  of  rocket  trajectory  when 
v'  » v0.  Hence,  as  long  as  the  flrat-stage  motion  covers  a ground  range  of  500  n.mi.  or  less, 
or  A'  < 10*,  we  have  a sufficient  condition  for  the  validity  of  the  flat  earth  approximation  for 
the  first  stage.  This  condition  is  very  often  satisfied  by  first-stage  motion,  although  in  our  equa- 
tions it  does  not  have  to  be  satisfied  by  higher  stages  of  a multistage  rocket. 


It  is  possible  to  go  beyond  the  flat  earth  approximation  for  the  first  stage,  and  we  already 
have  done  so  in  the  statement  of  Eq.  (25).  The  initial  problem  of  using  v„  - a»e  x r,  is  that 
we  don’t  know  ri  ahead  of  time,  since  it  is  determined  as  a subsequent  step  based  on  the  velo- 
city determination.  We  have  exploited  a flat  earth  approximation  to  substitute  for  ri  a vector 
from  the  center  of  the  earth  through  the  launch  point.  This  allowed  us  to  determine 
v - v„  + v‘,  where  r’  was  computed  without  the  effects  of  a rotating  earth.  This  itself  is  only 
an  approximation  for  the  contribution  Av*  i to  v',  although  we  shall  see  later  that  it  is  quite  a 
good  one.  With  v(/)  determined,  r(r)  followed  from  Eqs.  (27)  through  (34).  Next,  one  may 
choose  to  investigate  the  earth's  curvature  effect  for  the  first  stage  by  substituting  rt  from  this 
determination  into  x r,  as  the  start  of  another  iteration.  The  vector  rt  will  be 

specified  by  the  inertial  coordinates  yx,  Lu  and  /t.  From  then  on,  the  computations  of  Eqs. 
(26)  through  (34)  are  repeated,  with  some  modifications,  in  this  new  iteration.  The 
modifications  are  necessitated  by  the  new  orientation  of  the  vector  v„.  The  true  flight  path,  at 
least  the  one  determined  from  the  calculation  which  preceded  the  new  iteration,  is  shown  in 
Fig.  7 along  with  other  data  which  will  be  used  to  determine  the  aforementioned  modifications. 
Also  shown  is  the  flight  path  for  v'.  The  easterly  direction  at  (Lt,  /,)  is  the  same  as  that  at 
(£.  i.  /)).  We  will  use  the  given  information  about  the  spherical  triangle  ABC  in  Fig.  7 
(namely,  the  two  angles  lt  - /„  and  a0  and  included  side  ir/2  - Lfl)  to  determine  azimuth  angle 
a , and  ^i,  which  is  the  angular  distance  in  the  plane  of  r’  associated  with  the  first  stage.  These 
calculations  will  specify  the  orientation  of  r„.  From  standard  spherical  trigonometry  formulae 
[1],  we  have 


cos(/|  - /„)  + sin(/,  - /„)  cot  a„ 


or 


tan  L[  - (cos  Z,0)-1  Jsin  cos (lx  - (,)  + cot  a0  sin(/t  - (,)|. 
which  determines  £,j;  then 

sin  - sin(/t  - /„)  cos  Li /sin 


(35) 

(36) 


1! 


’■■-.■mw-ifPF''7 


i W*4l  «■•  tjt.  krtbuMbi 
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t ig  1 — Flight  path  projeiilmi  tin  nonrolatton  laiinulelongiiuite  gritl.  The 
diagram  make*  poaalblo  the  Inclu.tlon  of  earth's  rotation  tout  earth's  curvature 
fur  the  Ural  stage, 


determines  Azimuth  angle  a,  is  found  from  Ed.  (31): 

sin  at  - sin  a„  cos  /*/ cos  L\.  (37) 

Now  the  modiflestions  of  Eqs.  (27)  through  (34)  consist  of  the  following: 

val  - vu  cos  a i and  v„„  - v„  sin  a,  - (38) 

which  replaces  Eq.  (27)  - 

V*  - v‘  cos  + v0n  cos(*  - *|). 

V,  - v1  sin  t-  v(„t  sln($  - ^), 

and 

v - Ivj  + vi  + v|l),/,t  (39) 

which  is  used  instead  of  Eq.  (29);  the  other  equations  are  unchanged.  The  new  iteration  will 
result  in  new  velocities  and  positions,  which  may  be  made  the  basis  of  another  iteration,  if  it  is 
necessary  to  satisfy  convergence  criteria, 


dy  n/dt  w dy  f! di  — 2 o>  ^ x v^  — x (w  p x r), 

where  the  subscripts  /?  and  F refer  to  rotating  and  inertial  frames,  respectively.  When 
expressed  in  rotating  frame  coordinates,  the  aerodynamic  forces  associated  with  the  first  term 
on  the  right  properly  depend  only  on  rocket  motion.  The  latter  two  terms  are  Coriolis  and  cen- 


trifugal terms.  We  estimate  in  Appendix  C how  large  the  velocity  contributions  from  the  vari- 
ous forces  might  typically  be  for  a large  rocket  vehicle  at  cutoff  of  the  first  stage  (1-3].  In  order 
of  decreasing  values,  we  find  Avr:  Avc:  v0:  Av0:  AvCo,:  Avc<„  — 1000:  300:  ISO:  50:  10:  1.5, 
where  we  list,  in  order,  effects  of  thrust,  gravity,  initial  velocity  from  earth’s  rotation,  drag 
(which  is  approximately  cancelled  by  thrust  increase  due  to  ambient  pressure  decrease  [2,3]), 
Coriolis  force,  and  centrifugal  force.  These  estimates ‘indicate  the  conclusion  that  the  principal 
effect  of  a rotating  earth  is  to  add  in  the  velocity  v0  as  prescribed  in  Eqs.  (25)  and  (26)  and  dis- 
cussed subsequent  to  these  equations.  The  Coriolis  and  centrifugal  contributions  appear  to  be 
very  small,  so  that  Ava1  could  be  approximated  quite  well  by  the  approach  of  Section  2,  which 
disregards  the  effects  of  a rotating  earth  and  atmosphere. 

The  preceding  results  notwithstanding,  one  may  wish  to  improve  the  specification  of 
Avk1,  particularly  if  computational  simplicity  can  be  preserved.  The  accuracy  can  be  improved 
if  we  include  only  the  component  of  velocity  effect  from  Coriolis  and  centrifugal  forces  in  the 
particular  plane  of  motion  which  is  associated  with  the  disregard  of  rotational  effects.  To  illus- 
trate the  mathematical  point  that  this  is  a good  approximation  for  the  small  correction  from 
Coriolis  and  centrifugal  terms,  let  us  consider  what  happens  in  Eqs.  (29)  and  (30)  when  V = 0 
and  v„  « v'.  In  Eq.  (29)  we  have 

v - [v'2  + 2v'voU  cos  + v2J1/2  =v'  1 4-  cos  + 0 

and  in  Eq.  (29)  sin  i/f  = v'sin  In  these  equations  only  the  component  voM  enters  in  first 
order  for  v and  \ Jr.  The  only  first  order  correction  involving  v0l  that  can  be  found  is  in  the 
determination  of  a in  Eq.  (32),  and  this  is  not  important.  Hence,  along  with  Ruppe  [3],  we 
transform  to  the  rotating  frame  coordinates  in  Eqs.  (2)  and  (3)  by  setting 

t " 4>R  + to’ 

where 

to  ~ £e  cos  K sin  ao. 

When  Eq  (40)  is  substituted  into  Eqs.  (2)  and  (3),  and  the  steps  which  led  to  Sq3.  (6)  and  (7) 
are  redone,  it  is  found  that 

v*  - - sin  - — + ~ cos(a  + 0)  (41) 

Vkt  - - ~ rtl  ~ ~ cos  ip  + -jj  sin  (a  + 0)  + + 2vH  d>0,  (42) 

where 

yK  = (Av*j)f. 
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The  symbol  (Av*))f  has  been  defined  in  connection  with  Eq.  (26).  One  now  sees  additional 
Coriolis  and  centrifugal  terms  in  the  equation  of  motion.  In  the  same  way  that  Eqs.  (II)  and 
(12)  were  obtained,  one  finds 


v* 


sin  * + 


T 

M 


cos  a 


-U _ - a 


cos*  + ~|l  + 


bCt 

ba 


Spy  h 

IT 


sin  o + 2v**0 


(43) 

(44) 


In  terms  of  a programmed  deflection  function  *(r),  one  obtains  the  angle  of  attack  a 

from 


A/Vfi(»  - 2*J  + MUC/r1)  [1  - vjr/K  - r3»2/XI  cos  * 
T {1  + (8Ct/8a)  Spv*/2D 


(45) 


The  solution  for  vN  is  found  as  in  Eq.  (19).  It  is 

X ' ■ v*/v,  - " m (1  - t\)  - (\/ne)  (/,)t>  (46) 

where  ( and  A are  given  by  Eqs.  (15)  and  (17),  respectively,  and 

* - T/M0g'(R)  - r/^o  . (47) 

Here, 

g'(r)  - K/r1  - - g(r)  - r*02  (48) 

is  an  'apparent*  acceleration  of  gravity  and 

('.)<  - »i»  * <'{'  «» 

As  indicated  previously,  the  centrifugal  correction  to  gravity  is  only  a few  tenths  of  a percent. 
For  the  first  stage  xo  * 0 and,  from  Eq.  (39),  we  should  include  the  initial  condition 

*(0)  - 2*o  (first  stage)  (50) 

in  our  programmed  deflection  function.  This  is  similar  to  the  conclusion  reached  by  Ruppe  [3]. 
Since  - 7.292  x 10“5  rad/s,  it  is  seen  that  *(0)  in  Eq.  (50)  is  of  the  order  of  5 x 10~3 
deg/s,  which  corresponds  to  1*  in  200  s.  Similarly,  R*o  Is  on  the  order  of  10“ 2 m/s2.  These 
are  evidently  very  small  corrections,  in  accordance  with  our  previous  estimates,  so  that  AvX)  in 
Eq.  (24)  and  V in  Eq.  (26)  could  be  computed  to  a good  first  approximation  as  if  there  were  no 
earth's  rotation  effects. 

The  computational  procedure  for  relatively  simple  analyses  of  powered  rocket  ascent  is 
now  clear,  it  may  be  summarised  as  follows.  The  most  significant  part  of  the  rocket  trajectory 
is  assumed  to  be  above  the  atmosphere  and  associated  with  second  and  higher  stages,  with 
v0  « v\  so  that  for  any  qualitative  analysis  one  may  ignore  earth's  rotation  effects.  For  this 
purpose,  one  uses:  (a)  a programmed  deflection  function  *(/)  [2,3);  (b)  a thrust  given  by  Eq. 
(13)  with  correctly  chosen  exhaust  velocity  v,  [3,6),  (c)  Eqs.  (15)  through  (22)  for  the  veto* 
city,  altitude,  and  ground  range  of  the  rocket  in  ail  stages,  where  r - /?  is  an  expedient  approxi- 
mation for  low  altitude  rocket  flights  and  analytic  results  [3];  (d)  Eq.  (23)  for  checks  on  the 
smallness  of  the  angle  of  attack;  and  (e)  Eq.  (31)  for  the  determination  of  earth's  curvature 
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effects  in  the  computation  of  latitude,  azimuth  angle,  and  longitude  of  the  rocket  flight  path 
projection  on  the  surface  of  a nonrotating  earth  (hence,  an  inertial  latitude-longitude  grid).  If 
the  conditions  warrant  a flat  earth  approximation,  as  determined  by  the  ground  range  calcula- 
tion (e.g.,  if  4'  < 10^,  then  step  (e)  can  be  dispensed  with.  Also,  ii'  the  approximation  r =»  R 
has  been  made  in  step  (c),  one  may  improve  these  results  by  starting  with  them  in  another 
iteration  of  step  (c)  without  the  approximation. 

If  one  wishes  to  improve  quantitative  accuracy  by  taking  earth's  rotation  effects  into 
account,  he  should  And  the  following  summary  useful.  One  could  start  with  the  results  in  the 
preceding  paragraph  and  simply  add  in  v0  to  the  previous  velocities  v as  in  Eq's.  (25)  and  (26), 
where  rt  was  determined  in  the  previous  results.  Actually,  to  be  consistent,  one  should 
redetermine  the  flrst-stage  results  of  the  preceding  paragraph  with  the  use  of  Eqs.  (45)  through 
(50).  There  probably  will  be  no  appreciable  differences  from  the  previous  flrst-stage  results 
obtained  without  regard  for  earth's  rotation  effects,  but  the  results  should  at  least  be  checked. 
If  a flat  earth  approximation  is  valid  for  the  first  stage,  one  computes  component  and  total 
speeds  from  Eqs.  (27)  and  (29)  with  the  approximation  in  Eq.  (25),  heading  ip  from  Eq.  (30). 
and  azimuth  angle  a from  Eq.  (32).  The  primed  entities  in  these  equations  are  known  at  this 
point,  having  already  been  computed.  The  altitude  and  flight  path  projection  are  computed, 
presumably  through  numerical  integration,  in  Eq.  (34).  If  a flat  earth  aproximation  is  not  valid 
for  the  first  stage,  one  computes  component  and  total  speeds  from  Eqs.  (38)  and  (39),  where 
a l is  substituted  for  0|  (primed  entities  are  known).  Heading,  azimuth  angle,  altitude,  and 
flight  path  projection  are  again  determined  from  Eqs.  (30),  (32),  and  (34).  These  results  can 
then  be  used  as  the  start  of  another  iteration  for  obtaining  greater  accuracy  in  the  inclusion  of 
earth's  rotation  effects  when  the  flat  earth  approximation  for  the  first  stage  to  not  sufficiently 
accurate.  The  description  of  this  iteration  procedure  is  included  in  the  paragraph  which 
includes  Eq«.  (35)  through  (39).  It  is  expected  that  no  more  than  one  iteration  would  be 
required.  Even  one  iteration  might  not  be  worthwhile  in  view  of  the  basic  approximations  we 
have  made  regarding  thrust  and  aerodynamic  effects.  The  computations  are  greatly  simplified  if 
a flat  earth  approximation  can  be  made  for  the  total  powered  rocket  trajectory.  All  that  would 
be  required  for  a complete  specification  of  orbit  injection  conditions  (described  in  the  next  sec- 
tion) is  the  altitude  calculation  from  the  preceding  paragraph,  with  possible  alterations  which 
might  arise  from  Eqs.  (46)  through  (50)  for  the  first  stage  and  simple  vector  addition  of  veloci- 
ties, as  in  Eqs.  (25)  and  (26).  One  could  start  the  other  way  around,  however,  with  desired 
orbit  injection  conditions  (v,  \fi,  and  a)  and  obtain  the  launch  profile  conditions  (v\  0',  o^) 
needed  to  obtain  them;  here  v'  and  4>'  are  rocket  burnout  parameters.  This  method  is  explained 
in  Appendix  D. 

Finally,  if  one  wants  complete  quantitative  accuracy  regarding  all  the  effects  we  have 
mentioned— or  if  one  has  sufficient  information  about  aerodynamic  effects  and  other  factors  of 
the  system,  such  as  drag  and  lift  coefficients,  a model  atmosphere,  thrust  variability,  etc.,  and 
wants  to  include  these  effects  In  his  analysis— he  may  achieve  this  in  a full  numerical  solution 
to  the  problem.  The  discussion  of  Kooy  16]  is  particularly  helpful  in  this  regard.  Rocket  velo- 
city and  coordinates  can  be  straightforwardly  obtained  by  Kooy’s  procedure  (6)  in  a geocentric 
equatorial  coordinate  system  by  a Runge-Kutta  integration  method.  The  conversion  of  these 
results  to  latitude,  longitude,  azimuth,  altitude,  and  heading  variations  is  simple,  and  is  dis- 
cussed in  the  next  section.  Kooy’s  method  appears  to  be  quite  feasible,  as  well  as  accurate, 
with  the  use  of  modern  computers. 
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4.  ORBIT  TRAJECTORY  FROM  INJECTION 
CONDITIONS  AT  ROCKET  BURNOUT 

At  the  end  of  powered  rocket  ascent  (i.e.,  "burnout")  the  payload  is  assumed  to  enter  an 
elliptical  orbit  above  the  earth’s  atmosphere  in  the  spherical  gravitational  field  of  the  earth. 
This  satellite  motion  thus  satisfies  (in  an  inertial  frame) 

r + (A/r*)  r - 0,  (51) 

where  A - 1.407645  x 10u  =*  6.275  x 104  - is  the  same  parameter  as  used  pre- 
sec' sec' 

viously  (see,  e g.,  Eq.  (2)).  The  small  corrections  in  the  equation  of  motion  due  to  the  oblate 
earth,  aerodynamic  drag,  etc.,  can  be  handled  by  the  techniques  of  perturbation  theory  (1),  but 
will  not  be  included  in  the  present  report.  Dotting  Eq.  (51)  with  r and  integrating  with  respect 
to  time,  one  obtains  an  energy  constant  of  the  motion 

E - v2/2  - K/r.  (52) 


Similarly,  vector  multiplying  Eq.  (51)  by  r and  integrating  with  respect  to  time,  one  obtains  an 
angular  momentum  constant  of  the  motion: 

h-rxr.  (53)  1 


Hence,  the  motion  is  planar.  Now  if  one  vector  multiplies  Eq.  (51)  by  h and  integrates  it  with 
respect  to  time,  one  finds  [1,S] 

r x h - K^~  + e . (54) 

where  e is  a vector  integration  constant.  If  we  define  9 as  the  angle  between  e and  r,  the  dot 
product  of  this  equation  with  r yields  the  orbit  equation  {1,3) 

— - -^-(1  + e cos  9),  where  e - «*/ 1 + ~r~  (55) 

r h1  V A' 

is  the  eccentricity  of  the  orbit.  Fig.  8a  shows  the  elliptical  orbit  circumscribed  about  the  earth 
centered  at  one  of  the  foci,  and  Fig.  8b  shows  an  associated  vector  diagram.  Also  shown  are 
several  parameters  used  in  discussing  the  ellipse  and  the  motion  of  the  satellite.  The  parame- 
ters rf,  r„,  /,  and  a are  the  perigee,  apogee,  latus  rectum,  and  semimajor  axis  lengths,  respec- 
tively; their  values  can  be  derived  from  Eq.  (S5).  The  angle  9 is  measured  from  the  perigee 
vector  (in  the  direction  of  e),  as  shown  in  Fig.  8b,  and  the  velocity  and  heading  of  the  satellite 
are  shown  in  a manner  consistent  with  their  previous  use  in  this  report.  From  the  relation 

^ (56) 

dt  d9 


and  Eq.  (55),  one  obtains  (1) 


e sin  0 - (rv2/A)  sin  0 cos  < \> 
e cos  9 - (rv2/A)  cos2  0-1. 


These  are  useful  relations  for  many  purposes,  but  we  use  them  here  for  an  unambiguous  deter- 
mination of  the  injection  angle  9,  and  eccentricity  e in  terms  of  the  injection  conditions  r,,  v(, 
and  0,  at  burnout.  The  motion  of  the  satellite  in  orbit  is  thus  determined  from  injection  condi- 
tions, Eq.  (55),  in  which 
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if 
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E — v,ty2  - K/r,  and  h - r,  v,  cos 


(58) 


and  the  expression  for  the  time  f,(0)  it  takes  a satellite  to  travel  within  one  revolution  from 
perigee  to  the  angle  0.  This  is  (1] 


1»/J 


*<*>  - T 


{2  tan 


-t 


l - r 


t'/j 


1 + e 


$ 

tan  - 


e (1  - sin  0 
l + e cos  O 


(59) 


Hence,  for  example,  the  time  it  takes  to  travel  from  B,  to  B is 

t(9,  -0)  - r,(0)  - r,(0,). 


(60) 


The  only  other  thing  we  need  to  know  is  the  orientation  of  the  ellipse  with  respect  to  the 
earth  in  an  inertial  coordinate  system.  This  can  be  specified  from  the  vector  directions  of  h and 
e,  which  are  determined  from  injection  conditions.  Hence 

h - r,  x y, , (61) 

and  from  Eq.  (52) 

xh-~T--p  llv'J  “ t\  r<  " (r-  v)  TJ  (62) 
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The  inertial  coordinate  system,  which  we  have  previously  mentioned,  is  the  geocentric- 
equatorial  system  (1,5)  in  which  a latitude-longitude  grid  is  superimposed  upon  the  surface  of 
the  earth,  but,  unlike  the  conventional  latitude-longitude  grid  which  rotates  with  the  earth,  this 
grid  is  fixed  in  space  and  coincides  with  the  conventional  grid  at  the  time  of  launch  of  the 
rocket  from  the  earth.  The  latitude  and  longitude  changes  in  Eqs.  (31)  and  (34)  are  under- 
stood to  be  relative  to  this  space-fixed  latitude-longitude  grid.  A picture  of  satellite  motion  in 
the  geocentric-equatorial  system  is  shown  in  Fit.  9.  The  Cartesian  coordinate  directions  l ),  k 
(unit  vectors)  form  a right-hsnded  set,  with  K pointing  north  (perpendicular  to  the  equatorial 
plane)  and  / pointing  in  the  vernal  equinox  direction  of  the  sun.  A vector  V has  Cartesian 
components  ( Vx . Fv,  Vt)  in  this  system  or  components  ( F,,  F#.  F4)  in  the  associated  spheri- 
cal coordinate  system.  This  vector  can  also  be  represented  in  circular  perifocal  (S)  coordinates 
associated  with  the  orbital  plane  (see  Fig.  8)  as  ( Vr,  Vt,  Ft),  where  i is  in  the  direction  of  h. 
The  transformation  between  the  latter  two  coordinate  iystems  is 

f - f (63) 

P - - cos  a & + sin  a ♦ 

£<-fxP«-sin«6-cosa£, 


where  a symbol  with  a caret  denotes  a unit  vector  in  the  direction  indicated.  Here  angle  a in 
Fig.  9 is  the  azi.nuth  angle  of  satellite  motion,  as  described  previously.  The  angle  B:  is  the  co- 
latitude of  satellite  motion;  i.e.,  in  torms  of  latitude  L 


« - w/2  - L . (64) 

wheie  L > 0 in  the  northern  hemisphere  and  L < 0 in  the  southern  hemisphere.  To  com- 
plete the  transformation  to  Cartesian  coordinate  axes  l ),  k , we  use  the  expression  for  the 
azimuthal  angle  ♦ 

* - (0)  + f,  (65) 

where  <t»,  (0)  is  the  azimuth  of  the  Greenwich  meridian  at  the  time  of  rocket  launch.  The  long- 
itude of  the  satellite  / is  positive  if  it  is  east  of  the  Greenwich  meridian,  negative  if  west.  Now 
from  simple  trigonometry. 
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f - cos  Z cos  ♦ / + cos  Z sin  d>  i + sin  Z.  £ 

6 - sin  Z cos  ♦ / + sin  Z sin  ♦ J - c o$  L ft 

4 - - sin  ♦ / + cos  ♦ ).  (66) 


This  completes  the  orthogonal  transformations  between  the  above  coordinate  systems.  Hence, 
for  exempt  from  Eq.  (63) 

f - (-  sin  a sin  Z cos  ♦ + cos  a sin  4]  / 

+ l-  sin  a sin  Z • - cos  a cos  #1  J (67) 

+ sin  a cos  L ft, 

which  can  be  calculated  from  the  orbit  injection  conditions  for  a,  Z,  and  L 


Also  shown  in  Fig.  9 are  other  angles  used  to  describe  the  orientation  of  the  ellipse  in  the 
geocentric  equatorial  system.  These  are  the  inclination  /,  the  longitude  of  the  ascending  node 
O,  and  the  argument  of  perigee  /9.  The  inclination  angle  is  the  angle  between  ft  and  h,  found 
from  Eqs.  (63)  and  (64)  to  be 

cos  i - i -ft  - - sin  a £ ft  - sin  a cos  Z.  (68) 

In  general,  0 < # < w radians,  but  this  equation  restricts  inclination  angles  to  the  range 
|Z|  < / < n - |Z|,  depending  on  a.  For  inclinations  0 < / < w/2,  we  have  generally  east* 
erly  motion,  or  so  called  "direct*  orbits  (S).  For  w/2  < i < w,  we  have  'retrograde*  orbits  (5). 
If  we  define  the  line  of  nodes  as  a line  from  the  center  of  the  earth  through  the  point  of  cross- 
ing of  the  satellite  through  the  equatorial  plane,  and  if  the  unit  vector  A has  this  direction,  then 

i£(h  S )/  sin  i 

— — : ( cos  Z cos  a f + sin  Z cos  a 0 - sin  Z sin  a ♦!.  (69) 

sin  / 

This  crn  be  evaluated  from  orbit  injection  conditions. 


(6?): 


Before  we  can  finish  the  specification  of  U r.rd  fi,  we  shall  need  to  calculate  e from  Eq. 


» - ~ v x h - .*  - ~ cos  4 - 1 


vA 
K 


f - sin  6 


vA  . , 

v * ■ i 


f + sin  # Ifos  a 8 - sin  a 4>). 

fV 


which  can  b«.  evaluated  from  orbit  injection  conditions.  For  fl  one  obtains 
cos  0 - A / - (sin  /)“'  / • (ft  x I)  - - (sin  <)_l  f • J 

— (sin  a sin  Z sin  ♦ + cos  a cos  d>] , 
sin  i 

where  0 < 0 < 2w  rad  in  general,  but 

n < w if  A • ) > 0, 


(70) 


(71a) 
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or  its  equivalent,  if 


Hence,  from  Eq.  (67) 


i 1 > 0. 


ft  < tr  tf  — sin  a sin  L cos  ♦ + cos  a sin  $ > 0. 
One  also  obtains  from  Eqs.  (69)  and  (70) 

cos  0 - ii  • e 


1 


e sin  i 


(cos  l.  cos  a 


~ cos  # - 1 1 + sin  L ~ sin  + , 


where  0 < 0 < 2w  in  general,  but  if  e • k > 0,  0 < w,  I.e., 


0 < ir  if  ~ cos  # - I 


v/| 

sin  Z.  — — sin  + cos  a cos  JL  >0. 


(7  lb) 


(72a) 


(72b) 


AH  the  entities  in  Eqs.  (68)  through  (72)  are  calculable  from  known  orbit  injection  conditions, 
and  the  specification  of  the  elliptical  orbit  is  complete.  For  simplicity,  one  might  choose  to  set 
4>f  (0)  - 0 and  d>  - / in  the  preceding  expressions,  which  amounts  to  a redefinition  of  the  / 
and  } directions. 


5.  VEHICLE  FLIGHT  SEEN  BY  INERTIAL 
AND  EARTH  OBSERVERS 

In  Sections  2 and  3 the  functions  for  the  inertial  frame.  Lit),  lit), 
yit),  ait),  vit),  and  4(r),  were  specified  for  rocket  flight.  In  Section  4 additional  informa- 
tion was  derived  to  specify  these  entities  for  the  orbital  flight.  We  shall  specify  the  first  three 
of  these  along  with  the  counterpart  functions  LK(t),  /*(r),  and  y*it)  which  describe  what  the 
earth  observer  sees. 


It  is  very  simple  to  deal  with  the  period  of  powered  rocket  ascent,  i.e.>  the  time  frame 
0 < t < r„  where  f,  is  the  time  after  launch  to  burnout  or  orbit  injection.  In  Section  3 we 
have  discussed  the  computation  of  altitude,  latitude,  and  longitude  in  an  inertial  frame,  i.e., 
>{r).  Lit),  and  lit),  respectively,  during  this  time  period;  the  solution  for  yg it)  and  Lgit)  is 
thus  also  found,  since 


yRit)~ yit).  Lgit) -Lit).  (73) 

The  function  lg  it)  is  simply  found  from 

/*(/)  •lit)  -»Et.  (74) 

where 

mt  - 7.292  x 10’s  rad/s  - 15.04  deg/h;  (75) 

see  Eq.  (35)  and  [1.5], 

For  the  time  frame  t > the  satellite  is  in  the  elliptical  orbit  trajectory  discussed  in  the 
preceding  section.  The  angle  B in  Fig.  9 is  increasing  in  the  direction  of  motion  of  the  satellite, 
by  definition,  so  a computational  procedure  is  to  pick  a set  of  increasing  B values,  starting  with 
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9,  determined  from  Eq,  (57)  end  orbit  injection  conditions,  snd  to  determine  sn  associated  set 
of  times  from  Eqs,  (59)  and  (60),  starting  with  / — r.  From  spherical  trigonometry  formulae 
on  the  right  spherical  triangle  ABC  in  Fig.  9,  the  taw  of  sines  determines  latitude  L(t)  - L*(f) 
at  these  times  from 

sin  L(r)  • sin  i sin  l/J  t-  •(/)!.  (76) 

where  #(r)  is  given,  and  0 and  t have  previously  been  computed  from  orbit  injtcton  conditions 
by  Eos.  (61)  and  (72)  A cotangent  formula  (I)  determines  ((.)  et  these  times  Dorn 

tan  (/(/)  > ♦,((?)  - (l!  - tan  (0  + •(#>)  cos  t,  (77) 

where  0 is  determined  from  injection  conditions  by  Eq.  (7)),  and  (f)  is  given  by  Eq.  (72). 
Finally,  altitude  y(f)  is  determined  from  the  known  orbit  in  Eq.  (S3)  by  the  formula 

V(/)  -*0)  -rl#(r)J  - *(i).  (71) 

where  R{t)  is  the  radius  of  the  earth,  which  can  be  taken  as  approximately  constant,  and 
r (•(()!  for  the  orbit  is  given  by  Eq.  (55). 

«.  DISCUSSION 

In  Sections  2 and  3 this  report  develops  the  equations  required  for  an  analysis  of  powered 
rocket  trajectory  Most  of  the  effort  has  been  concentrated  on  the  approximations  which  permit 
the  simpler  analysis  techniques;  the  approach  to  the  more  complex,  completely  quantitative 
numerical  methods,  however,  has  at  least  been  indicated.  The  computational  procedures  for 
various  levels  of  simplicity  and  approximation  are  delineated  by  summary  discussions  near  the 
end  of  Section  3.  The  investigator  can  thus  choose  the  procedure  which  (Its  his  needs.  It  was 
recognised  that  a radial  gravitational  field  introduces  centrifugal  terms  in  an  inertial  frame, 
terms  which  arc  omitted  from  some  literature  sources  or  given  incorrectly.  A systematic  pro- 
cedure for  taking  earth's  rotation  effects  (including  atmosphere  corotation)  into  account  was 
also  described,  a subject  often  treated  incompletely  in  the  literature.  Perhaps  the  most 
significant  contribution  of  this  report  is  the  treatment  of  earth's  curvature  effects  and  how  the 
earth's  rotation  effects  are  integrated  with  it  in  the  simpler  analyses  of  powered  rocket  trajec- 
tory. The  literature  treatments  of  approximate  analyses  of  powered  rocket  trajectory,  at  least 
those  seen  by  the  author,  stay  within  the  confines  of  the  flat-earth  approximation;  this  approxi- 
mation is  discussed  in  Section  3.  The  full  numerical  solution  of  the  rocket  problem  (6]  does 
typically  include  earth's  curvature  effects. 

One  of  the  major  approximations  of  simpler  analysis  techniques  of  powered  rocket  trajec- 
tory is  the  cancellation  of  drag  and  thrust  increase  effects.  This  approximation  appears  to  be 
quite  good  for  many  rockets  (2,31,  i.e.,  within  2%  for  strategic  choice  of  the  constant  thrust 
value  for  each  stage  (3).  For  large  rockets,  however,  the  thrust  increase  effect  is  expected  to 
more  than  counterbalance  the  drag  effect,  so  that  one  may  wish  to  include  these  effects  expli- 
citly for  greater  quantitative  accuracy.  The  simplification  from  constant  thrust  is  then  lost  from 
the  analysis,  but  it  is  eliminated  anyway  when  angles  of  attack  become  large,  as  they  frequently 
do  in  higher  stage  motion.  One  may  then  wish  to  solve  equations  (8)  and  (9)  in  Section  2,  or 
equations  (10)  and  (II),  by  numerical  integration  (e  g.,  by  the  Runge-Kutta  method).  Infor- 
mation about  thrust  increase  and  aerodynamic  effects  is  included  in  Appendix  E.  For  somewhat 
greater  understanding  and  control  of  the  launch  trajectory,  one  may  wish  to  assume  a profile 
a(r)  for  the  angle  of  attack  and  solve  for  the  flight  path  heading  # and  for  speed  in  the  course 
of  the  numerical  interaction,  a job  suited  for  a computer.  Previously,  we  had  suggested,  for 
simplicity,  an  assumed  6(f)  profile  in  the  calculations. 
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Section  4 relates  the  rocket  burnout  parameters  obtained  in  Sections  2 and  3 to  the 
specification  of  the  elliptical  orbit  subsequently  traversed  by  the  payload  and  to  the  point  on  this 
orbit  where  it  is  injected.  The  equations  are  developed  somewhat  more  clearly  and  completely 
than  usual  here,  and  it  is  hoped  that  the  orbit  determination  part  of  the  trajectory  analysis  is 
thus  facilitated. 

In  Section  S the  trajectory  coordinates  of  the  powered  rocket  ascent  and  orbital  flight  are 
determined  as  the  altitude  and  the  latitude  and  longitude  of  the  flight  projection  on  the  earth's 
surface  as  a function  of  time  after  launch.  Two  cases  are  considered:  (1)  the  trajectory  u seen 
by  a heavenly,  inertial  observer  on  a space-fixed  latitude-longitude  grid  on  the  earth,  and  (2) 
the  more  important  case  of  the  common  earth-fixed  observer  who  rotates  underneath  the  orbit 
and  sees  the  rocket  and  payload  from  a different  perspective.  Actually,  the  first  case  is  just 
obtained  as  an  intermediate  step  to  the  determination  of  the  second  case. 

There  are  other  perturbations  which  have  not  been  considered  in  this  report,  such  as  the 
effect  of  the  oblate  earth,  variations  in  aerodynamic  forces,  solar  radiation  pressure,  etc.  (1] 
Some  of  these  effects  alter  the  long  term  orbital  motion  of  satellites,  which  is  of  no  particular 
consequence  for  many  simple  analyses  or  for  this  report,  which  concerns  itself  with  powered 
rocket  ascent  and  short-term  orbit  motion.  Variable  winds  on  the  rocket  are  assumed  to  be 
corrected  for  in  flight  as  part  of  the  vehicle  steering  control. 
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Appendix  A 

LAGRANGIAN  METHOD,  INCLUDING  THE  PRESENCE  OF 
NONCONSERVATIVE  FORCES 

A system  of  particles  is  effectively  in  equilibrium  under  the  influence  of  all  applied  and 
inertial  forces.  The  applied  forces  may  be  further  classified  as  conservative  (derivable  as  the 
netative  gradient  of  a scalar  potential)  or  nonconservative.  We  may  define  the  forces  as  fol- 
lows: 

fit/, 

F,  — — - — : Conservative  force  on  Ah  particle 

F':  Nonconservative  force  on  Ah  particle 

Inertial  force  on  Ah  particle.  (Al) 

Newton’s  equation  of  motion  is 

F,  + F,  - m,x,  - 0.  (A 2) 

here  restated  as  particle  equilibrium  under  all  three  kinds  of  forces.  Since  each  particle  is  in 
effective  equilibrium,  an  arbitrary,  inAnitesimal,  virtual  displacement  of  all  the  particles  will 
involve  no  work.  This  is  known  as  d’Alembert’s  principle  of  virtual  displacements.  Hence 

£( F,  + Pi  - mi*,)  • 8x,  - 0.  (A3) 

I 

Because  of  constraints  in  the  problem  (e.g.,  boundaries,  rigid  body  constraints  of  the  particles, 
etc.),  particle  motion  will  depend  on  a smaller  set  of  coordinates  (e.g.,  angles  of  rotation,  the 
center  of  mass  coordinates,  etc.),  which  are  referred  to  [4]  as  independent  generalized  coordi- 
nates g|,  ....  qf.  Hence, 

Xj  - *,  («!.  *2 «/).  (A4) 

8x,  " H-  (A5) 

J 

Therefore, 

i - ii<«  lo  i/)- 

One  obtains 

X m'ij  • 8 x,  - £ • — lij 


(A6) 
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ft 


8ft 


- I 4 Zt  **’  - £ It  «tf| ««; 

; (<«  8ft  , 2 «ft  / 2 


It  is  also  true  that 


/ / > oft 

Hence,  if  we  include  the  kinetic  energies  and  conservative  forces  in  the  Lagrangian  L as 


i(  •••  i ft>  ft*  •••)  " 51 


qj,  .,.) 


where  i,  is  determined  by  differentiating  Eq.  (A4),  then  Eq.  (A3)  can  be  written  as 

d BL  BL 


I 

J 


. . . - F 

dt  Bqj  Bqj  ") 


Bq}  - 0, 


where 


(A7) 

(A8) 


Fa-  If; 


bx, 


(A9) 


defines  the  generalized  nonconservative  force  associated  with  the  coordinate  ft.  Since  the  ft 
coordinates  are  independent  and  the  displacements  8 ft  are  arbitrary,  each  of  the  square  brackets 
in  Eq.  (A8)  vanishes,  and  we  therefore  obtain  Lagrange’s  equations  suitably  generalized  to 
include  nonconservative  forces. 


Appendix  B 

DERIVATION  OF  THE  LAGRANGIAN  AND  GENERALIZED  FORCES 

FOR  THE  ROCKET  PROBLEM 

The  coordinate  variables  which  can  vary  independently  in  the  problem  are  the  coordinates 
( r , <f»)  of  the  center  of  mass  and  the  angle  of  attack  a.  Now 

x,  - r +£,,  (Bl) 

where  r is  the  radius  vector  to  the  center  of  mass  and  jo  / is  the  vector  displacement  from  the 
center  of  mass  to  the  i th  particle  mass  element  in  the  rocket.  The  kinetic  energy  Tis  given  by 

T - y £ m,{i  +£,)2  - y Mi 2 + £ y »»/p2  + | £ «£/|  * r, 
but  since  £ «,£,  — 0 by  definition  of  the  center  of  mass,  the  last  term  vanishes.  Conse- 

i 

quently,  we  have  the  familiar  separation  of  center  of  mass  motion  from  internal  motion.  Now 

£ \ 'W/P2  “ £ T J2  “ £ T «/*'  “*  W/  *jm * ««<P/)„ 

/ L i L i l 

“ £ 4'  ml(&k m *ln  ~ “*  fa/)/  (p/)» 

I L 

" T £ «*  p//  Pit  - «*  “/  Pil  P/J 

- -5-  at ' • I • at',  where  I - £ m,[p,2J^  --££,]. 

In  this  case  of  planar  motion,  where  at  — az  (z direction  normal  to  plane),  / "*  £ «/p2. 


Hence, 


r - y M(r2  + r2*2)  + y /a2. 


The  only  conservative  force  in  the  problem  is  gravity,  which  can  be  related  to  a potential: 

. 3/A 

,2  “ r2  r ~ r * 

rl  r I r 

Hence,  the  appropriate  Lagrangian  which  includes  conservative  forces,  is 
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From  Fig.  2: 

K - L(L,  + D,  + T()  • (r  +£()  - (L  + D + T)  • t 

i “r 

- L cos  - D sin  ^ + T sin(t|i  + a + 0).  (B4) 

Similarly, 

F;  - £( L + D,  + T,)  • - - (L  + D + T)  r • * 

m —L  sin  — D cos  + + T cos(tf»  + a + 0)  r.  (B5) 

Finally, 

^-I(L(  + K/+T()  ^(r+£,) 

#£,  — 2 sin  8a/2  (f  X£((a)l  - 8a[£  xj^^a)] 

K - I(L,  + K,  + T,)  • (£  x£(}. 

Now  we  use  the  information  that  lift  and  drag  act  as  if  concentrated  at  the  center  of  pressure 
and  thrust  acts  as  if  localized  at  the  center  of  combustion: 

F„  - (L  + D)  • (£  x |j)  + T • (£  x|,)  - /2U  cos  a + D sin  a)  - /,  T sin  0.  (B6) 
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Appendix  C 


VELOCITY  INCREMENT  ESTIMATES  FOR  VARIOUS  EFFECTS 
IN  FIRST-STAGE  ROCKET  MOTION 


The  ideal  velocity  (i.e.,  from  thrust  only)  of  a large  rocket  at  first-stage  cutoff  [1-3]  could 
be  estimated  (e.g.,  see  Eq.  (8)  of  text)  as 


Avr-  / ' <T/M)<it  m (T/A#)/,. 


(Cl) 


where  t\  is  the  burn  duration  of  the  first  stage,  and  M is  an  average  mass.  The  gravity  loss 
from  this  velocity  is 

Av c =»*/,=»  0.3 (T/M)tx,  (C2) 

where  we  indicate  a 30%  gravity  loss  [2].  Now  Avr  might  be  =»  3050  m/s,  so  that  Avc  =*  900 
m/s.  The  drag  toss  from  A vr  is  probably  =»  5%  overall,  although  the  drag  force  can  reach  a 
maximum  of  about  25%  of  the  thrust  force  during  the  first-stage  motion  [2]. 


Hence, 


Av0  s*  150  m/s  s*  0.05  Av, , 


(C3) 


which  is  almost  canceled  by  thrust  increase  due  to  ambient  pressure  decrease  in  the  ascent 
[2,3].  By  comparison  we  have  a relatively  [arge  value  [5]  of  velocity  imparted  by  earth’s  rota- 
tion 


(C4) 


v0  *■  UfR  — 450  m/s  =*  0.15  Avr. 

We  estimate  the  Coriolis  efTect  as 

A Vfo,  =*  1 2«£V  dt  an  2a>g(T/M)  f^'ldt 

or 

Avfw  a»  wfr,  Avr  =*  0.01  Avr  (C5) 

for  a bum  time  t\  somewhat  more  than  2 min.  The  centrifugal  effect  contribution  is  given  as 


Av*,,,  =»  w|/lr,  at  <*st\y0  * 0.0015  Av7 


(C6) 
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Appendix  D 

DETERMINATION  OF  LAUNCH  PROFILE  PARAMETERS 
FROM  DESIRED  ORBIT  INJECTION  CONDITIONS 
IN  A FLAT  EARTH  APPROXIMATION 

If  we  are  given  a set  of  desired  orbit  iivjection  conditions,  such  as  altitude  y,  rocket  speed 
v,  heading  0,  sod  azimuth  angle  a,  we  face  the  question  of  the  launch  profile  cond"  as  needed 
to  obtain  them.  We  design  the  rocket  to  attain  the  speed  v'  and  y'  at  burnout  ^’’regarding 
earth's  rotation)  for  a final  heading  0',  and  need  to  determine  the  apparent  launch  azimuth  a0. 
The  relationship  between  these  parameters  is  obtained  with  the  help  of  the  vector  diagram  in 
Fig.  Dl.  The  notation  used  here  has  been  defined  in  Section  3 of  the  text  (cf.  also  Figs.  4 and 
6).  The  following  results  may  readily  be  obtained. 


N 


FI*.  Dl  — Vector  diagram  Tor  determination  of  launch  profile  parameters 
from  desired  orbit  injection  conditions  (flat  earth  approximation) 

v'  - I*  - V0|  - |v„  + v*  - v„|  - {v*  + (v,,  - v0)2}'/2 
- |vJ  - 2vv0  cos  0 sin  a + 

With  v'  thus  determined,  we  use  the  following  from  Fig.  Dl: 

sin  0’  - v sin  0/v\ 


(Dl) 


(D2) 


S7W£EP*Jjm.  v>- 


riJ^7^r*7"7?^?£^  '^T' 
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This  determines  •!/',  and  now  if  we  define 

Aa  - a - a0  (D3) 

and  use  the  law  of  sines  in  Fig.  Dl, 

sin  Aa  - v#  cos  a/(v'  cos  ifr'),  (D4) 

from  which  a„  is  determined.  Since  v0  is  in  the  horizontal  plane,  we  have 

/ - y.  (D4) 


A similar  set  of  equations  has  been  derived  by  Ruppe  [3]. 


Appendix  E 

THRUST  INCREASE  AND  AERODYNAMIC  EFFECTS 

For  a rocket  with  swivel  control  motors  the  total  thrust  includes  pressure  forces  and  is 
normally  presented  (2,3,6]  as 

T - Mvf  + A[pt  - p(y)].  (El) 

where  M is  the  mass  toss  rate  of  exhaust  gases,  v,  is  the  actual  average  axial  speed  of  these 
gases  relative  to  the  rocket,  A is  the  exit  aperture  area,  p,  is  the  pressure  of  the  exhaust  gases 
at  the  exit  aperture,  and  p(y)  is  the  ambient  pressure  of  the  atmosphere  at  altitude  y.  Nor- 
mally, p(y)  varies  from  p(o)  to  zero  for  first-stage  motion,  so  that  the  thrust  varies  from  its 
sea  level  value  Tsi  to  its  vacuum  value  Twc.  Hence,  T,\  and  7^  are  usually  given  for  the  first 
stage,  and  TMC  is  given  for  the  higher  stages  of  a rocket.  In  effect  the  ratio 

x s fw/r„  (E2) 

is  given  for  the  first-stage  motion,  and  it  is  easily  shown  that  (El)  can  be  rewritten  as 

T - 7^  [l  - — -44  • (E3) 

x p(o) 

Very  frequently  xis  the  range  1.12  to  1.15,  and  the  increase  of  thrust  with  altitude  is  thus  exhi- 
bited. 

As  we  have  indicated  in  the  text,  drag  tends  to  counterbalance  the  effect  of  thrust 
increase.  The  aerodynamic  forces  of  lift  L and  drag  D are  important  in  first-stage  motion, 
where  angles  of  attack  a are  normally  programmed  to  be  small  (e.g.,  a "=*"  3s).  From  [2]  one 
then  has  for  the  aerodynamic  forces: 

D - (Coo  + CDL  «21  (pvJ/2)  S 

L - [dQ/da]  a (pvJ/2)  S,  (E4) 

where  p is  the  atmospheric  density,  v is  the  rocket  speed,  and  5 is  an  appropriate  reference 
cross-sectional  area  (e.g.,  that  for  the  first  stage)  to  which  the  aerodynamic  coefficients  are 
referred.  These  coefficients  for  a two-stage  rocket  are  shown  by  Ehricke  (2]  in  his  Fig.  5-8  as  a 

function  of  local  mach  number,  which  is  the  speed  v divided  by  the  speed  of  sound  at  the  alti- 

tude of  the  rocket.  One  may  thus  infer  the  aerodynamic  forces  from  standard  atmosphere  data 
and  the  aerodynamic  coefficients. 
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